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Abstract. We review the various theories which have been proposed along the years to explain 
the origin of the stellar initial mass function. We pay particular attention to four models, namely 
the competitive accretion and the theories based respectively on stopped accretion, MHD shocks 
and turbulent dispersion. In each case, we derive the main assumptions and calculations that 
support each theory and stress their respective successes and failures or difficulties. 
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1. Introduction 

Stars are the building blocks of our universe and understanding their formation and 
evolution is one of the most important problems of astrophysics. Of particular importance 
is the problem of the initial mass function of stars (Salpeter 1955, Scalo 1986, Kroupa 
2002, Chabrier 2003) largely because the star properties, evolution and influence on the 
surrounding interstellar medium strongly depend on their masses. It is generally found 
that the number of stars per logarithmic bin of masses, dN/dlogM, can be described by 
a lognormal distribution below 1 Mq, peaking at about ~0.3 Mq, and a power-law of 
slope —1.3 for masses between 1 and 10 (e.g. Chabrier 2003). It should be stressed 
that the IMF of more massive stars is extremely poorly known. 

Several theories have been proposed to explain the origin of the IMF, invoking vari- 
ous physical processes that we tentatively classify in four categories: theories based on 
recursive fragmentation or pure gravity (e.g. Larson 1973), theories based on pure statis- 
tical argument, invoking the central limit theorem (e.g. Zinnecker 1984, Elmegreen 1997, 
Adams & Fatuzzo 1996), theories based on accretion and, finally, theories invoking the 
initial Jeans mass in a fluctuating environment. We will focus on the two latter ones, 
which appear to be favored in the modern context of star formation. 

2. Theories based on accretion 

2.1. Competitive accretion 

The theory of competitive accretion has been originally proposed by Zinnecker (1982) 
and Bonnell et al. (2001). It has then been used to interpret the series of numerical 
simulations similar to the ones performed by Bate et al. (2003). 

The underlying main idea is that the accretion onto the star is directly linked to its 

f It should be kept in mind that the initial mass function is measured only up to about 8 
Mq, in young stellar clusters, and is inferred only indirectly for larger masses (see e.g. Kroupa 
2002) 
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mass in such a way that massive stars tend to accrete more efficiently and thus become 
disproportionally more massive than the low mass stars. The accretion rate is written as: 

M, = TipVrelRlcc. (2-1) 

where p is the gas density, Vrei is the relative velocity between the star and the gas while 
Race is the accretion radius. Bonnell et al. (2001) consider two situations, namely the 
cases where the gravitational potential is dominated by the gas or by the stars. 

2.1.1. Gas dominated potential 

Let R be the spherical radius, p the gas density and the number density of stars. 
In the gas dominated potential case, Bonnell et al. (2001) assume, following Shu (1977), 
that the gas density profile is proportional to R~^ . They further assume that rt, oc R^^ . 
The accretion radius. Race is assumed to be equal to the tidal radius given by 

/ M \^^^ 

Rudai ^ 0.5 — ^ R, (2.2) 

\ -^^^enc J 

where Mcnc is the mass enclosed within radius R. This choice is motivated by the fact 
that a fluid particle located at a distance from a star smaller than Rudai is more sensitive 
to the star than to the cluster potential and will thus be accreted onto this star. 

The mass of gas within a radius R, M{R), is proportional to R (since p (x R~^). The 
infall speed is about ~ y^GM{R) / R, and, assuming that the stars are virialized, one 
gets Vrei — Vin- The number of stars, diV,, located between R and R + dR is given by the 
relation dN^, — n,(i?) x AnR^dR cx dR. Thus can 12.11 combined with the expression of 
Vrei and Race Icads to the relation M* c>c (Af,/i?)^/'^ and, after integration, oc R~^, 

— 1/2 

R oc M, . Consequently, we obtain: 

dN (X IVC^^^dM^. (2.3) 

Even though the mass spectrum is too shallow compared to the fiducial IMF, dN/dM oc 
M~'^-^, it is interesting to see that this power law behaviour can be obtained from such 
a simple model. Note, however, that the model implies that stars of a given mass are all 
located in the same radius, which points to a difficulty of the model. 

2.1.2. Star dominated potential 

When the potential is dominated by the stars located in the centre of the cloud, the 
density is given by p oc i?^"^/^, which corresponds to the expected density distribution 
after the rarefaction wave has propagated away (Shu 1977). The velocity is still assumed 
to be Vrei oc R~^/^. The accretion radius is now supposed to correspond to the Bondi- 
Hoyle radius as the gas and the star velocities are no longer correlated. This leads to 
Af* = npVreiRBH whcrc Rbh oc AI^,/V^^i. It follows 

M* oc A/~2. (2.4) 

One can then show (Zinnecker 1982, Bonnell et al. 2001) that under reasonable assump- 
tions, dN oc M^dM^,. This estimate is in better agreement with the Salpeter exponent, 
although still slightly too shallow. These trends seem to be confirmed by the simulations 
performed by Bonnell et al. (2001) which consists in distributing 100 sink particles in a 
cloud of total mass about 10 times the total mass of the sinks initially (their Fig. 3). 

2.1.3. Difficulties of the competitive accretion scenario 

As obvious from the previous analytical derivations, finding an explanation for the 
Salpeter exponent with the competitive accretion scenario appears to be difficult, even 
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though numerical simulations (as the ones presented in Bonnell et al. 2001) seem to suc- 
cessfully achieve this task. However, although the IMF exponent is close to the Salpeter 
one in the star dominated potential case, this scenario entails by construction the Bondi- 
Hoyle accretion, which is at least a factor 3 lower than the mass infall rate resulting 
from gravitational collapse at the class O and I stages and leads to too long accretion 
times compared with observations (Andre et al. 2007, 2009). Another difficulty of this 
scenario is that it does not explain the peak of the IMF which might be related to the 
Jeans mass (see §[3]). Finally, it is not clear that competitive accretion can work in the 
case of non-clustered star formation, for which the gas density is much too small. As no 
evidence for substantial IMF variation among different regions has yet been reported, 
this constitutes a difficulty for this model as a general model for star formation. Perhaps 
this scenario applies well to the formation of massive stars. 



The principle of this type of models is to assume that the accretion of gas onto the stars 
or the dense cores is a non-steady process, stopped because of either the finite reservoir 
of mass or the influence of an outflow which sweeps up the remaining gas within the 
vicinity of the accreting protostar. 

The flrst studies were performed by Silk (1995) and Adams & Fatuzzo (1996). They 
flrst relate the mass of the stars to the physical parameters of the cloud such as sound 
speed and rotation and then assume that an outflow whose properties are related to the 
accretion luminosity stops the cloud collapse. Using the Larson (1981) relations, they can 
link all these parameters to the clump masses. Since the mass spectrum of these latter 
is known (e.g. Heithausen et al. 1998), they infer the IMF. 

A statistical approach has been carried out by Basu & Jones (2004). These authors 
assume that the dense core distribution is initially lognormal, justifying it by the large 
number of processes that control their formation (and invoking the central limit theorem). 
Then, they argue that the cores grow by accretion and postulate that the accretion 
rate is simply proportional to their mass, M = 7M — M{t) = Moex.p{'-ft), leading to 
log Af = fi = fj,o +jt. Finally, they assume that accretion is lasting over a flnite period of 
time given by f(t) = S exp(—St). The star mass distribution is thus obtained by summing 
over the accretion time distribution. 



where a = S/j and (Tq characterizes the width of the initial dense core distribution. As S 
and a are controlled by the same types of processes, their ratio is expected to be of the 
order of unity and thus f{M) exhibits a powerlaw behaviour close to the fiducial IMF. 

In a recent study, Myers (2009) develops similar ideas in more details, taking into ac- 
count the accretion coming from the surrounding background. Adjusting two parameters, 
he reproduces quite nicely the observed IMF (his figure 5). 

A related model has also been developed by Bate & Bonnell (2005) based on an idea 
proposed by Price & Podsiadlowski (1995). They consider objects that form by fragmen- 
tation within a small cluster and are ejected by gravitational interaction with the other 
fragments, which stops the accretion process. Assuming a lognormal accretion rate and 
an exponential probability of being ejected, these authors construct a mass distribution 
that can fit the IMF for some choices of parameters. 

In summary, the stopped accretion scenario presents interesting ideas and, providing 



2.2. Stopped accretion 




(2.5) 
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(typically 2 or more) adequate adjustable parameters, can reproduce reasonably well the 
IMF. However, the very presence of such parameters, which characterizes our inability to 
precisely determine the processes that halt accretion, illustrates the obvious difficulties 
of this class of models, and their lack of predictive power and accuracy. 



3. Gravo-turbulent theories 

While in the accretion models, turbulence is not determinant, it is one of the essential 
physical processes for the two theories presented in this section, although the role it 
plays differs in both models, as shown below. The theories proposed along this line 
seemingly identify cores or pre-cores and are motivated by the strong similarity between 
the observed CMF and the IMF (e.g. Andre et al. 2010). 

The first theory which combined turbulence and gravity was proposed by Padoan et 
al. (1997). In this paper, the authors consider a lognormal density distribution - density 
PDF computed from numerical simulations (e.g. Vazquez- Semadeni 1994, Kritsuk et al. 
2007, Schmidt et al. 2009, Federrath et al. 2010) are indeed nearly lognormal - and select 
the regions of the flow which are Jeans unstable. By doing so, they get too stiff an IMF 
(typically dN/dM cx M^^) but nevertheless find a lognormal behaviour at small masses, 
a direct consequence of the lognormal density distribution, and a powerlaw one at large 
masses. 



3.1. Formation of cores by MHD shocks 

The idea developed by Padoan & Nordlund (2002) is slightly different. These authors 
consider a compressed layer formed by ram pressure in a weakly magnetized medium. 
They assume that the magnetic field is parallel to the layer and thus perpendicular to 
the incoming velocity field. The postshock density, pi, the thickness of the layer. A, 
and the postshock magnetic field, Bi, can be related to the Alfvenic Mach numbers, 
■Ma = v/va {v is the velocity and Va the Alfven speed), and preshocked quantities, po 
and Bq according to the shock conditions: 

pi/po ^ Ma , A/L ~ , Bi/Bo - Ma , (3.1) 

where L is the scale of the turbulent fiuctuation. Note that for classical hydrodynamical 
isothermal shocks, the jump condition is typically c>c M"^. The dependence on Ma instead 
oi M^ stems from the magnetic pressure which is quadratic in _B. As we will see, this is 
a central assumption of this model. 

The typical mass of this perturbation is expected to be 

m ~ piA^ ~ poXa (^-^) ^PoL^M-^ (3.2) 

As the flow is turbulent, the velocity distribution depends on the scale and u ~ i", 
with a — {n — 3)/2, E(k) cx fc~" being the velocity powerspectrunlj. Combining these 
expressions with eqn p.2|) . they infer 

m^^(^y'\ (3.3) 



Ma,0 \Lo 

where ip is the largest or integral scale of the system and Ma.o the corresponding Mach 
number. To get a mass spectrum, it is further assumed that the number of cores, N{L), 



f n is denoted /3 in Padoan & Nordlund (2002), more precisely n — 2 — P 
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formed by a velocity fluctuation of scale L, is proportional to i Combining this last 
relation with eqn p.3p leads to 



For a value of n — 3.74 (close to what is inferred from 3D numerical simulations), one 
gets N{ni) ~ m~^-^^, very close to the Salpeter exponent. 

So far, gravity has not been playing any role in this derivation and the mass spectrum 
that is inferred is valid for arbitrarily small masses. In a second step, these authors 
consider a distribution of Jeans masses within the clumps induced by turbulence. As 
the density in turbulent flows presents a lognormal distribution, they assume that this 
implies a lognormal distribution of Jeans lengths and they multiply the mass spectrum 
(|3.4p by a distribution of Jeans masses, which leads to 



The shape of the mass spectrum stated by eqn 13.51 is very similar to the observed IMF 
(see for example the figure 1 of Padoan &; Nordlund 2002). 

Note, however, that difhculties with this theory have been pointed out by McKee 
& Ostriker (2007) and Hennebelle & Chabrier (2008). Eqn p.l|) . in particular, implies 
that in the densest regions where dense cores form, the magnetic field is proportional 
to the density, in strong contrast with what is observed both in simulations (Padoan & 
Nordlund 1999, Hennebelle et al. 2008) and in observations (e.g. Troland & Heiles 1986). 
This is a consequence of the assumption that the magnetic field and the velocity field are 
perpendicular, which again is not the trend observed in numerical simulations. In both 
cases, it is found that at densities lower than about 10'^ cm"-^, B depends only weakly 
on n while at higher densities, B cx ^/p. This constitutes a problem for this theory, 
as the index of the power law slope is a direct consequence of eqn (13. ip . Assuming a 
different relation between B and p, as the aforementioned observed one, would lead to 
a slope stiffer than the Salpeter value. Furthermore, the Salpeter IMF is recovered in 
various purely hydrodynamical simulations (e.g. Bate et al. 2003), while the Padoan & 
Nordlund theory predicts a stiffer distribution (dN/dM oc AI~^) in the hydrodynamical 
case. Another important shortcoming of this theory is that it predicts that turbulence, 
by producing overdense, gravitationally unstable areas, always promotes star formation, 
while it is well established from numerical simulations that the net effect of turbulence 
is to reduce the star formation efficiency (e.g. MacLow & Klessen 2004). 



Recently, Hennebelle & Chabrier (2008, 2009, HC08, HC09) proposed a different theory 
which consists in counting the mass of the fluid regions within which gravity dominates 
over the sum of all supports, thermal, turbulent and magnetic, according to the Virial 
condition. In this approach, the role of turbulence is dual: on one hand it promotes star 
formation by locally compressing the gas but on the other hand, it also quenches star 
formation because of the turbulent dispersion of the flow, which is taken into account in 
the selection of the pieces of fluid that collapse. 

The theory is formulated by deriving an extension of the Press & Schechter (1974, 
PS) statistical formalism, developed in cosmology. The two major differences are (i) the 
underlying density field, characterized by small and Gaussian fluctuations in the cosmo- 
logical case while lognormal in the star formation case, and (ii) the selection criterion, 
a simple scale- free density threshold in cosmology while scale-dependent, based on the 
Virial theorem in the second case. That is, fluid particles which satisfy the criterion (see 



N{m)d\ogm ~ m 



3/(s-")dlogm. 



(3.4) 




(3.5) 



3.2. Turbulent dispersion 
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HC08) 

(y,L>+3(CfO'<-i?pot/Af (3.6) 

are assumed to collapse and form a prestellar bound core. The turbulent rms velocity 
obeys a power-law correlation with the size of the region, the observed so-called Larson 

relation, (F^s) = x ( if^j > with Vq ~ Ikms"^ and -q ~ 0.4-0.5 (Larson 1981). 

The principle of the method is the following. First, the density field is smoothed at a 
scale i?, using a window function. Then, the total mass contained in areas which, at scale 
i?, have a density contrast larger than the specified density criterion (5^, is obtained by 
integrating accordingly the density PDF. This mass, on the other hand, is also equal to 
the total mass located in structures of mass larger than a scale dependent critical mass 
M^, which will end up forming structures of mass smaller than or equal to this critical 
mass for collapse (see HC08 §5.1). 

pe^]i{5)VB{5)d5 ^ / M' N{M') P{R,M')dM' . (3.7) 



In this expression, M'j^ is the mass which at scale R is gravitationally unstable, = 
log{p'j^/p) and = M^/ (CmR^), Cm being a dimensionless coefficient of order unity. 
Vr is the (turbulent) density PDF, assumed to be lognormal, while P{R, M') is the 
conditional probability to find a gravitationally unstable mass, M' embedded into Mf^ at 
scale R, assumed to be equal to 1 (see HC08 §5.1.2 and App. D). Note that this expression 
is explicitly solving the cloud in cloud problem as the mass which is unstable at scale 
R is spread over the structures of masses smaller than (right hand side integral). 
Therefore, by construction, all the gravitationally unstable regions in the parent clump 
that will eventually collapse to form individual prestellar cores are properly accounted 
for in this theory (see HC08 §5.1). 

Taking the derivative of eqn (j3.7p with respect to R, we obtain the mass spectrum 

While the second term is important to explain the mass spectrum of unbound clumps 
defined by a uniform density threshold (as the CO clumps), it plays a minor role for 
(Virial defined) bound cores and can generally be dropped (see HC08 for details). 
After soma algebra and proper normalisation, one gets 

1 1 + (1 - J])MlR^'i m\ ' ^ exp(-(TV8), 



N{M) = 2Afo ^ ^ " * 1 X X ^ (3-9) 

R^ [I + {2r] + l)MlR^'i] \Ry v^<J 

where R ^ R/\°, 3]^ = In |(1 + MlR'^'^) / li^^ , JVq = p/M° and M% denote the 

usual thermal Jeans mass and Jeans length, respectively, and 

M{R) = M/M°j ^R{1+ MlR^"^) (3.10) 

denotes the unstable mass at scale R in the turbulent medium. The theory is controlled 
by two Mach numbers, namely 

x-^4^f^yMo.-i,o,f^)7^)-, (3.„) 



VsCs Vlpc/ VO-lpc/ V0-2kms 
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defined as the non-thermal velocity to sound speed ratio at the mean Jeans scale Aj 
(and not at the local Jeans length), and the usual Mach number, A4, which represents 
the same quantity at the scale of the turbulence injection scale, Li, assumed to be the 

/y2> 1/2 

characteristic size of the system, M. = ^-^ — . 

The global Mach number, M, broadens the density PDF, as ct^ = ln(l + b^Ai^), 
illustrating the trend of supersonic turbulence to promote star formation by creating 
new overdense collapsing seeds. 

The effect described by A^, is the additional non thermal support induced by the tur- 
bulent dispersion. In particular, at large scales the net effect of turbulence is to stabilize 
pieces of fluid that would be gravitationally unstable if only the thermal support was 
considered. This is illustrated by eqn p.lOp which reduces to M = i? when A^* = 0. 
In particular, for a finite cloud size, the gas whose associated turbulent Jeans length is 
larger than the cloud size is not going to collapse. 

When AJ* <^ 1, i.e. the turbulent support is small compared to the thermal one, 
eqn p.9|) shows that the CMF at large masses is identical to the Padoan et al. (1997) 
result, i.e. dN/d\og M cx M~^. On the other hand, when Mt, ~ 1, dN/d\og M cx 
^-(Ti+i)/(2n-4)^ where the index of the velocity powerspectrum n is related to 77 by the 
relation -q = (n — 3)/2 (see HC08). As n ~ 3.8 — 3.9 in supersonic turbulence simulations 
(e.g. Kritsuk et al. 2007), turbulent dispersion leads to the correct Salpeter slope and 
Larson velocity-size relation. 

Comparisons with the Chabrier (2003) IMF have been performed for a series of cloud 
parameters (density, size, velocity dispersion) and good agreement has been found (IIC09) 
for clouds typically 3 to 5 times denser than the mean density inferred from Larson (1981) 
density-size relation. Comparisons with numerical simulations have also been performed. 
In particular, Schmidt et al. (2010), performing supersonic isothermal simulations with 
various forcing, have computed the mass spectrum of cores supported either by pure 
thermal support or by turbulent plus thermal support. Their converged simulations show 
very good quantitative agreement with the present theory, confirming that turbulent 
support is needed to yield the Salpeter index. Note that Schmidt et al. (2010) use for 
the density PDF the one they measure in their simulations which is nearly, but not 
exactly lognormal. Comparisons with the results of SPH simulations (Jappsen et al. 
2005) including self-gravity and thermal properties of the gas have also been found to be 
quite successful (HC09). 

3.3. Difficulties of the gravo-turbulent theories 

One natural question about any IMF theory is to which extent it varies with physical 
conditions. Indeed, there is strong observational support for a nearly invariant form 
and peak location of the IMF in various environments under Milky Way like conditions 
(see e.g. Bastian et al. 2010). Jeans length based theories could have difficulty with the 
universality of the peak position, since it is linked to the Jeans mass which varies with the 
gas density. Various propositions have been made to alleviate this problem. Elmegreen 
et al. (2008) and Bate (2009) propose that the gas temperature may indeed increase 
with density, resulting in a Jeans mass which weakly depends on the density, while HC08 
propose that for clumps following Larson relations, there is a compensation between the 
density dependence of the Jeans mass and the Mach number dependence of the density 
PDF, resulting in a peak position that is insensitive to the clump size. 

A related problem is the fact that massive stars are often observed to be located in 
the densest regions, where the Jeans mass is smaller. Indeed, Mj cx p~^/^ when a purely 
thermal support is considered, whereas Mj cx p^^ when turbulence is taken into account 
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(assuming that V cx L°-^) (see HC09). This constitutes a difficulty for theories based on 
Jeans mass although, as seen above, the issue is much less severe when turbulent support 
is considered as massive stars can be formed at densities only few times smaller than the 
densities at which low mass stars form. Another possibility is that dynamical interactions 
between young protostars may lead to the migration of massive stars in the center of the 
gravitational well. 

Furthermore, the dependence of the freefall time on the Jeans mass should also modify 
the link between the CMF and the IMF, as pointed out by Clark et al. (2007). This 
is particularly true for theories which invoke only thermal support. When turbulent 
support is included, the free-fall time is found to depend only weakly on the mass, with 
tff (X M^/"* (see McKee & Tan 2003 and HC09 App. C), resolving this collapsing time 
problem. We stress that this time represents the time needed for the whole turbulent 
Jeans mass to be accreted. It is certainly true that, within this turbulent Jeans mass, 
small structures induced by turbulent compression will form rapidly. Their total mass, 
however, is expected to represent only a fraction of the total turbulent Jeans mass because 
the net effect of turbulence is to decrease the star formation efficiency (see IIC08). 

Generally speaking, the fragmentation that occurs during the collapse could consti- 
tute a problem for theories invoking Jeans masses. Although this problem is far from 
being settled, it should be stressed that such a fragmentation process is not incompatible 
with the calculations performed by HC08. As shown by eqn 13. 7[ the presence of small 
self-gravitational condensations induced by turbulence at the early stages of star forma- 
tion and embedded into larger ones is self-consistently taken into account in the theory. 
Moreover, the SPH simulations performed by Smith et al. (2008) show a clear correlation 
between the initial masses within the gravitational well and the final sink masses up to 
a few local freefall times (see Chabrier & Hennebelle 2010 for a quantitative analysis), 
suggesting that the initial prestellar cores do not fragment into many objects. As time 
goes on, the correlation becomes weaker but seems to persist up to the end of their run. 
Massive stars, on the other hand, are weakly correlated with the mass of the potential 
well in which they form. Whether their mass was contained into a larger more massive 
well with which the final sink mass would be well correlated remains an open issue, which 
needs to be further investigated. At last, both the magnetic field (e.g. Machida et al. 2005, 
Hennebelle & Teyssier 2008) and the radiative feedback (Bate 2009, Offner et al. 2009) 
will reduce the fragmentation, suggesting that the core-sink correlation found in Smith et 
al. (2008) should improve if such processes were included. Clearly, these questions require 
careful investigations. 

4. Conclusion 

We have reviewed the most recent theories which have been proposed to explain the 
origin of the IMF. Due to limited space, it was not possible to cover all of them (e.g. Kunz 
& Mouschovias 2009). Two main categories received particular attention: the theories 
based on accretion and the ones based explicitly on turbulence. It should be stressed that 
these theories are not all exclusive from each other and may apply in different ranges 
of mass. For instance, the turbulent dispersion theory calculates the distribution of the 
initial mass accretion reservoirs; it is not incompatible with the stopped accretion theories 
and with the competitive accretion as long as mass redistribution/competition occurs 
within one parent core reservoir. The question as to whether one of these mechanisms is 
dominant is yet unsettled. Detailed comparisons between systematic sets of simulations, 
as done in Schmidt et al. (2010), or observations, and the various analytical predictions 
is clearly mandatory to make further progress. 
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